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Abstract 

The present paper provides a primer on confirmatory factor 
analysis. An overview of the basic concepts and processes of 
confirmatory factor analysis will be presented. Finally, 
Campbell, Gillaspy, and Thompson's (in press) confirmatory factor 
analysis of the Bern Sex-Role Inventory (BSRI) (Bern, 1981) will be 
used to illustrate the application of CFA in the behavioral 
sciences. 
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A Primer on Confirmatory Factor Analysis 
Conventional factor analysis can be defined as a set of 
analytic techniques, “designed to examine the covariance 
structure of a set of variables and to provide an explanation of 
the relationships among those variables in terms of a smaller 
number of unobserved latent variables called factors” (Daniel, 
1988, p. 2) . As such, factor analysis provides the researcher a 
means of reducing large numbers of variables to a more manageable 
size and of investigating the latent constructs underlying the 
variables. Ker linger (1986) explained that: 

Factor analysis serves the cause of scientific 
parsimony. It reduces the multiplicity of tests and 
measures to greater simplicity. It tells us, in 
effect, what tests or measures belong together ... It 
thus reduces the number of variables with which the 
scientist must cope. (p. 569) 

Factor analytic methods have long been used to examine the 
validity of psychological constructs and the psychometric 
instruments used to measure those constructs. Nunnally (1978, 
p. Ill) noted that construct validity has often been referred to 
as factorial validity.” Thus, Nunnally (1978, pp. 112-113) 
argued that “factor analysis is intimately involved with 
questions of validity,” and that “Factor analysis is at the heart 
of the measurement of psychological constructs. This proposition 
is echoed by Kerlinger (1986, p. 468) , who stated that, "Factor 
analysis is perhaps the most powerful method of construct 
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validation." Thompson and Daniel (in press) review the history 
of applications of factor analysis in validity studies, and 
recommend selected best practices. 

Exploratory Factor Analysis (EFA) 

There are two major types of factor analytic techniques; 
conventional or exploratory factor analysis (EFA) and 
confirmatory factor analysis (CFA) . The distinction between 
exploratory and confirmatory factor analysis is based on the 
purpose of the data analysis. As implied by the name, 
exploratory factor analysis is used to explore the data when the 
researcher does not have sufficient evidence to form a hypothesis 
about the number or nature of the factors underlying the data. 
Confirmatory factor analysis, on the other hand, is a technique 
in which the number and the composition of factors are specified 
prior to the analysis or extraction of factors. Historically, 
the majority of factor analytic studies have been exploratory 
(Gorsuch, 1983; Kim & Mueller, 1978). Exploratory factor 
analytic methods have primarily been used in the development of 
theory or in the definition of theoretical constructs. 

Unfortunately, the use of exploratory factor analysis as a 
means of evaluating the construct validity of test scores or of 
serving the purposes of scientific parsimony is inappropriate in 
some circumstances. Mulaik (1988) argued that the extraction of 
factors does not provide "justified or authenticated knowledge" 
(p. 267). 

There are several features of exploratory factor analysis 
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that inherently limit its usefulness for evaluating construct 
validity. First, the factor structures yielded in exploratory 
factor analysis are determined by the mechanics of the method. 

For example, different extraction procedures, such as Cattell’s 
scree test (Cattell, 1966) or Guttman’s eigenvalue greater than 
one rule (Guttman, 1954) , can produce different numbers of 
factors (Zwick & Velioer, 1986) . 

Even in exploratory methods, researchers impose certain 
assumptions on the data that do not always honor the 
relationships between the variables. For example, orthogonal 
rotations require that the factors be perfectly uncorrelated, 
while oblique rotations require that the factors be correlated 
with one another. In other words, when using either orthogonal 
or oblique rotation methods, the researcher a priori may impose 
potentially mistaken assumptions onto the data. Thus, it should 
be clear that the factors obtained during an exploratory factor 
analysis are very much dependent on specific theories and 
mechanics of extraction and rotation procedures. Despite 
misconceptions to the contrary, exploratory factor analysis does 
not simply allow the data to “speak for themselves.” 

How to interpret the factors once they have been extracted 
during EFA is another problem. Nunnally (1978) noted that 
interpretation of factors measured by few variables is often 
complicated. Mulaik (1972) also raised concern about 
interpreting the results of exploratory factor analysis; 

This difficulty most often comes about because the 
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researcher lacks even tentative prior knowledge 
about the processes which produce covariation among 
the variables studied and has no basis on which to 
make his interpretations. In these circumstances 
the interpretations given the factors may be 
nothing more than tautological transformations of 
the names of the original variables. Difficulty is 
also encountered when the factors obtained 
represent confounded effects and the researcher is 
unable to decide which of these effects is unique 
to the factor — a problem which may come about from 
random selection of variables. (p. 363) 

Due to these and other shortcomings, Mulaik (1987, p. 301) 
cautions that "It is we who create meanings for things in 
deciding how they are to be used. Thus we should see the folly 
of supposing that exploratory factor analysis will teach us what 
intelligence is, or what personality is." 

Confirmatory Factor Analysis fCFA) 

Confirmatory factor analysis offers the researcher a more 
viable method for evaluating construct validity. In confirmatory 
factor analysis, a model positing the number and the composition 
of the factors is determined prior to the analysis and tested 
against the data in hand. Thus, confirmatory factor analysis 
enables the researcher to explicitly test hypotheses concerning 
the factor structure of the data. For this reason, Gorsuch 
(1983, p. 134) stated that, "Confirmatory factor analysis is the 
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more theoretically important-and should be the much more widely 
used“Of the two major factor analytic approaches." 

Basic Procedures of Conf irmatory Factor Analysis 

Confirmatory factor analyses can be conducted using computer 
programs such as LISREL VII (Joreskog & Sorbom, 1989) . The first 
step in a confirmatory factor analysis is to propose competing 
models or hypotheses about the factor structure of variables 
based on theory or existing data. Different models of the data 
are proposed by "fixing” and "freeing" specific parameters, i.e., 
the factor coefficients, the factor correlation coefficients, and 
the variance /covariance of the error of measurement. Fixing a 
parameter refers to setting the parameter at a specific value 
based on one's expectations. Thus, in fixing a parameter the 
researcher does not allow that parameter to be estimated in the 
analysis. Within the command computer file, such as the one 
presented in Appendix A, fixed variables are indicated by being 
zero values. 

Freeing a parameter refers to allowing the parameter to be 
estimated during the analysis by fitting the model to the data 
according to some theory about the data. The competing models or 
hypotheses about the structure of the data are then tested 
against one another. 

The confirmatory factor analysis yields several different 
statistics for determining how well the competing models 
explained the covariation among the variables or fit the data. 
These statistics are referred to as "fit statistics.” These 
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indicators of model fit include the chi square /degrees of freedom 
ratio, the Bentler Comparative Fit Index (CFI) (Bentler, 1990) , 
the parsimony ratio, and the Goodness-of-f it Index (GFI) 

(Joreskog & Sorbom, 1989) . Using the chi square/degrees of 
freedom ratio, a model is considered to have a better fit with 
the data when the ratio is 2-to-l or lower. Using the GFI, 
better models are those in which the GFI approaches 1.00. 

The parsimony ratio is also important when interpreting the 
data. This ratio takes into consideration the number of 
parameters estimated in the model. The fewer number of 
parameters needed to specify the model, the simpler or more 
parsimonious is the model. In science, parsimony is sought to 
explain phenomena because more parsimonious solutions are more 
likely to be true and thus typically more generalizable. The 
parsimony ratio can also be multiplied by a fit statistic to 
produce an index of both the overall efficacy of the model to 
explain the covariance among the variables and the parsimony of 
the proposed model. With the variety of fit statistics available 
to the researcher for evaluating model fit, Campbell, Gillaspy, 
and Thompson (1995, p. 6) advocated that "fit should be 
simultaneously evaluated from the perspective of multiple fit 
statistics. ” 

An Example of Conf irmatorv Factor Analysis 

To illustrate the basic concepts and procedures of 
confirmatory factor analysis, Campbell, Gillaspy, and Thompson's 
(in press) study of the covariance structure of the scores of the 
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Bern Sex-Role Inventory (BSRI) (Bern, 1981) is reviewed. The Bern 
Sex-Role Inventory is a popular measure of androgyny. The 
development of the BSRI can be traced to Constantinople (1973) 
who argued that persons could possess stereotypically masculine 
and stereotypically feminine psychological traits in any 
combination, regardless of physical gender. For example, persons 
who are both masculine and feminine in their psychological 
outlook are termed "androgynous". In the BSRI, masculinity and 
femininity are considered to be two separate constructs. Prior 
to Constantinople's work, masculinity and femininity were 
considered to be opposite poles of a single construct. 

The BSRI yields scores on a masculine (M) and feminine (F) 
scale. These scores are then dichotomized, creating a 2X2 
contingency table with four cells: (a) masculine (high M, low F) ; 
(b) feminine (high F, low M) ; (c) androgynous (high M, high F) ; 
and (d) undifferentiated (low M, low F) . Figure 1 presents the 
BSRI 2X2 contingency table. This scoring implies a model in 
which the two constructs or factors, masculine and feminine, are 
uncorrelated . 



INSERT FIGURE 1 ABOUT HERE 

The structure underlying responses to the Bern Sex-Role 
Inventory has been investigated using various analytic methods 
across diverse samples (see Thompson, 1989) . Thompson and 
Melancon (1986) provide an example of the use of exploratory 
methods with scores from the measure. Confirmatory methods have 
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been applied to BSRI data from adolescents (Thompson & Melancon, 
1988) . Second-order confirmatory methods have also been used 
(Marsh, 1985) . 

The present study investigated the covariance structure 
associated with the BSRI data provided by 791 graduate and 
undergraduate students enrolled at a large university. The 
sample was predominately white (82.9%), though the sample also 
included Hispanics (9.5%), and African-Americans (4.2%). There 
were slightly more women (50.9%) in the sample. The mean age was 
20.23 (SD=4.04). 

Three completing models of the factors underlying the BSRI 
were evaluated: 

Model #la (v=40) . This model posited a single bipolar 
factor defined by the 40 variables or items of the 
BSRI. This model tested the plausibility of the view 
of masculinity and femininity prior to Constantinople 
(1973) . The expectation was that 20 items measuring M 
(or F) should have positive factor parameters, while 
the remaining 20 items should have negative parameters. 

Model #2a (v=40) . This model posited two uncorrelated 
factors defined by the 40 variables — 20 items per 
factor. This was the model implied by conventional 
scoring of the BSRI that creates the four-fold 
classification typology. 

Model #3a (y=40) . This model posited two correlated factors 
defined oy the 40 variables — 20 items per factor. This 
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model presumed that there are two discernable M and F 
dimensions, but either (a) that the dimensions may have 
a stable but non-zero correlation or (b) that the 
correlation of the dimensions may vary over samples or 
time, even though the M and F constructs themselves 
remain invariant. 

As noted by Thompson (19f;4) the characteristics of 
reliability and validity inure to scores and not to tests. In 
fact, sometimes scores from shorter tests are more reliable than 
scores from longer tests (Thompson, 1990, p. 586) . The scores on 
the BSRI provide a powerful example of this situation. 
Specifically, the 20-item short-form of the BSRI generally yields 
comparable or more reliable scores (a^ ranging from .84 to .86; 

Up ranging from .84 to .87) than does the 40-item long-form (a^ 
ranging from .86 to .87; Cp ranging from .75 to .78), especially 
on the Feminine scale (Bern, 1981, p. 14) . Therefore, parallel 
models of the structure underlying only the responses to the 
subset of the 20 BSRI short-form were also evaluated: 

Model #lb (y=20) . This model posited a single bipolar 
factor defined by the 20 variables. 

Model #2b (Y=20) . This model posited two uncorrelated 
factors defined by the 20 variables — 10 items per 
factor. 

Model #3b (y=20) . This model posited two correlated factors 
defined by the 20 variables — 10 items per factor. 

Each model was specified by freeing (a) one factor parameter 
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per variable, (b) the factor correlation coefficients, and (c) 
the measurement error variance for each variable. All other 
parameters were fixed. Table 1 presents a comparison of 
exploratory vs confirmatory models as regards the parameters that 
are estimated. 

INSERT TABLE 1 ABOUT HERE 

Confirmatory factor analyses were conducted using LISREL 
covariance structure analyses (Joreskog & Sorbom, 1989) . 

Bivariate correlation matrices were used as the basis for each 
analysis, to produce "scale-free" parameters. Correlation 
matrices could be used because all the models involved variables 
correlating with only one factor, and each factor had factor 
variance fixed to one (Cudeck, 1989) . For most other CFA 
applications, the variance/covariance matrix must be analyzed, as 
Cudeck (1989) explained. 

Table 2 presents the fit statistics (Bentler, 1990, 1994) 
for the six models. Table 3 presents the maximum- likelihood 
parameter estimates for Model #3b. 

INSERT TABLES 2 AND 3 ABOUT HERE 

Several features of this analysis are noteworthy. First, 
neither a model positing no factors nor a model positing a single 
bipolar factor fit the data, as indicated by the various fit 
statistics reported in Table 2. The failure to fit a bipolar 
single factor supports Constantinople's (1973) original theory as 
regards these constructs. 
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Second, the correlation (r = -.076) between the two factors 
in Model #3b (v=20) was negligible, as reported in Table 3. The 
correlation between the two factors in Model #3a (v=40) was -022. 
These results suggest that the two constructs may be essentially 
orthogonal, at least in the present study, as implied by a 
classification scheme presented as the 2X2 contingency table 
typically employed by researchers using the Bern Sex-Role 
Inventory . 

Third, the fit statistics presented in Table 2 would not 
make one completely sanguine about the validity of scores from 
the long form of the measure, or at least of scores computed 
using conventional scoring keys. Models #2a and #3a fit the data 
equally well, but neither provided a particularly good fit. 

The Table 2 fit statistics also suggest the superiority of 
scores from the short form. Again, researchers can not simply 
assume that scores from long forms inherently have better 
measurement characteristics than scores from shorter test 
versions. Consistent with previous literature (cf. Bern, 1981), 
the 20-item short-form of the Bern yielded comparable or more 
reliable scores (a^ = .82, Op = .89) for our data than did the 
40-item long-form (Om = .85, Cp = ‘SI), especially on the 
Feminine scale. 

Models #2b and #3b fit the data reasonably well. Therefore, 
given both the CFA and alpha coefficients, it is suggested that 
scores on the short-form of the BSRI may have more utility for 
the purposes of future research. 
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Suroroarv 

The present paper has attempted to demonstrate that 
confirmatory factor analysis is, in most cases, the tool of 
choice for evaluating construct validity. Confirmatory factor 
analysis holds a distinct advantage over exploratory factor 
analytic methods in that CFA enables the researcher to test 
competing hypotheses about the number and composition of the 
factors underlying the data. Exploratory methods, through the 
use of specific extraction and rotation methods, impose certain 
assumptions on the data that may not honor the nature and 
relationship among the variables. An example of a confirmatory 
factor analysis of the Bern Sex-Role Inventory (Bern, 1981) was 
used to illustrate the basic concepts and procedures of CFA. 
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Figure 1 

Bea Sex-Role Inventory 2X2 Contingency Table 



Masculine 



High 



High 



Feminine 

Low 



Androgynous 


Masculine 


Feminine 


Undifferentiated 
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Table 1 

Factor Matrices of (a) Exploratory Factor Analysis and (b) 
Confirmatory Factor Analysis for the Two Factor Hypothesis of the 
Short-form of the Bern Sex-Role Inventory 



(a) Exploratory Factor Analysis (b) Confirmatory Factor 

Analysis 



Items I 

Dimension 


II 


Items 


I 


II 




1 


Xii 


X12 


1 


Xii 


0 


M 


2 


^21 


X22 


2 


X21 


0 


M 


3 


^31 


^32 


3 


X31 


0 


M 


4 


X41 


X42 


4 


X41 


0 


M 


5 


^51 


X52 


5 


^51 


0 


M 


6 


^61 


X(S2 


6 


X(si 


0 


M 


7 


X71 


X?2 


7 


X71 


0 


M 


8 


^81 


^»2 


8 


^81 


0 


M 


9 


X91 


X92 


9 


X91 


0 


M 


10 


^101 


^;02 


10 


^101 


0 


M 


11 


^111 


^112 


11 


0 


^112 


F 


12 


^121 


X122 


12 


0 


X122 


F 


13 


^131 


^132 


13 


0 


Xl 32 


F 


14 


^141 


^142 


14 


0 


^142 


F 


15 


^151 


^152 


15 


0 


^152 


F 


16 


^161 


^162 


16 


0 


^162 


F 


17 


^171 


Xi 72 


17 


0 


^172 


F 


18 


^181 


^182 


18 


0 


^182 


F 


19 


^191 


Xi 92 


19 


0 


^192 


F 


20 


^201 


X202 


20 


0 


X202 


F 



X = the factor pattern/ structure coefficients; M; masculine 
items; F: feminine items. 
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Table 2 

Fit Statistics for Six Models 
(n = 791; y = 40 or 20) 



Model 



Statistic 


ilSL 


t2& 


ils. 


ilh 


#2b 


i2h 


v 


40 


40 


40 


20 


20 


20 


Null chi sq 


12628.22 


12628.22 


12628.22 


6093.77 


6093.77 


6093.77 


Null df 


780 


780 


780 


190 


190 


190 


Noncentral ity 


11848.22 


11848.22 


11848.22 


5903.77 


5903.77 


5903.77* 


Model chi sq 


8891.26 


5415.64 


5415.35 


2636.03 


926.40 


922.96 


Model df 


740 


740 


739 


170 


170 


169 


Noncentrality 


8151.26 


4675.64 


4676.35 


2466.03 


756.40 


753.96* 


NC / df 


11.01522 


6.31843 


6.32794 


14.50606 


4.44941 


4.46130*’ 


GFI 


0.482 


0.728 


0.728 


0.627 


0.884 


0.884 


Pars Ratio 


0.90244 


0.90244 


0.90122 


0.80952 


0.80952 


0.80476* 


GFI*Pars 


0.43498 


0.65698 


0.65609 


0.50757 


0.71562 


0.71141' 


CFI 


0.31203 


0.60537 


0.60531 


0.58230 


0.87188 


0.87229* 


Pars Ratio 


0.94872 


0.94872 


0.94744 


0.89474 


0.89474 


0.88947^ 


CFI*Pars 


0.29603 


0.57433 


0.57349 


0.52100 


0.78010 


0.77588* 



•Noncentrality = " df 

‘’Noncentrality / df 

‘Parsimony Ratio = Model df / [(variables * (variables + 1)) / 2] 
‘‘GFI * Parsimony Ratio 

crvT = r fNull x! - Null df) - (Model x! z Model df) 1 

(Null - Null df) 

^Parsimony Ratio = Model df / [(variables * (variables -I)) / 2] 
*CFI * Parsimony Ratio 



Note . These fit statistics are described by Bentler (1990, 1994) and by 
Mulaik, James, van Alstine, Bennett, Lind, and Stilwell (1989). Both il 
models ("a" and ”b”) posited a single bipolar dimension; the #2 models 
posited two uncorrelated factors; the #3 models posited two correlated 
factors. 
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Table 3 

Correlated Two-Factor Model 
Maximum Likelihood Parameter Estimates 



LAMBDA X 

DEFENDMY 

INDEPEND 

ASSERTIV 

STRONGPE 

FORCEFUL 

LEADERSH 

TAKERISK 

DOMINANT 

TAKESTAN 

AGRESSIV 

AFFECTIO 

SYMPATHE 

SENSITIV 

UNDERSTA 

COMPASSI 

SOOTHEHU 

WARM 

TENDER 

LOVECHIL 

GENTLE 

PHI 

MASCULIN 

FEMININE 



MASCULIN 

0.469 

0.394 

0.700 

0.640 

0.474 

0.591 

0.479 

0.663 

0.614 

0.628 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 



MASCULIN 

1.000 

-0.076 



FEMININE 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.651 

0.677 

0.684 

0.635 

0.771 

0.649 

0.770 

0.794 

0.383 

0.757 



FEMININE 

1.000 
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Appendix A 

Annotated LISREL Syntax File for BSRI Confirmatory Factory 

Analysis 



00168 


SUBTITLE 'A CREATE CORK MATRIX N=791 V*20' 






00169 1 


PRELIS 






00170 


/VARIABLES 






00171 


defendroy to agresslv affectio to gentle (CO) 






00172 


/TYPE-CORRELATION 






00173 


/MATRIX-OUT (COR) 






00174 ; 


LISREL 






00175 


/"O. null model **VCOV MATRIX** 


UNCOR 


FACTORS 


00176 


/DA NI«20 NO«791 MA=KM 






00177 


/MATRIX>IN(COR) 






00178 


/MO NX>20 NK-20 PH*ZE LX*SY TD»ZE 






00179 


/FR LX(1,1) LX(2,2) LX(3,3) LX(4,4) LX(5,5) 






00180 


/FR LX(6,6) LX(7,7) LX(8,8) LX(9,9) LX(10,10) 






00181 


/FR LX(11,11) LX(12,12) LX(13,13) LX(14,14) LX(15,15) 




00182 


/FR LX(16,16) LX(17,17) LX(18,18) LX(19,19) LX(20,20) 




00183 


/OU SE MR FS SL=1 TM=1200 ND=7 






00184 


LISREL 






00185 


/"I. FACTOR ANALYSIS BEM.DTA **CORR MATRIX** 


1 FACTOR" 


00186 


/DA NI«20 NO»791 MA=KM 






00187 


/MATRIX-IN (COR) 






00188 


/MO NX=20 NK=1 PH=SY TD=SY 






00189 


/LK 






00190 


/MASCULUN FEMININE 






00191 


/PA LX 






00192 


/I 






00193 


/I 






00194 


/I 






00195 


/I 






00196 


/I 






00197 


/I 






00198 


/I 






00199 


./I 






00200 


/I 






00201 


/I 






00202 


/I 






00203 


/I 






00204 


/I 






00205 


/I 






00206 


/I 






00207 


/I 






00208 


/I 






00209 


/I 






00210 


/I 






00211 


/I 






00212 


/PA PH 






00213 


/o 






00214 


/VA 1.0 PH(1,1) 






00215 


/OU SE MR MI FS SL«1 TM-1200 ND»3 






00216 


LISREL 






00217 


/"2. FACTOR ANALYSIS BEM.DTA **CORR MATRIX** 


UNCOR 


FACTORS 


00218 


/DA NI«20 NO«791 MA>KM 






00219 


/MATRIX«IN(COR) 






00220 


/MO NX«20 NK«2 PH»SY TD«SY 






00221 


/LK 






00222 


/MASCULUN FEMININE 






00223 


/PA LX 
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00224 /I 0 

00225 /I 0 

00226 /I 0 

00227 /I 0 

00228 /I 0 

00229 /I 0 

00230 /I 0 

00231 /I 0 

00232 /I 0 

00233 /I 0 

00234 /O 1 

00235 /O 1 

00236 /O 1 

00237 /O 1 

00238 /O 1 

00239 /O 1 

00240 /O 1 

00241 /O 1 

00242 /O 1 

00243 /O 1 

00244 /PA PH 

00245 /O 

00246 /O 0 



00247 /VA 1.0 PH(1,1) PH(2,2) 

00248 /OU SB MR HI FS SL=1 TH=1200 ND=3 

00249 LISREL 

00250 /"3. FACTOR ANALYSIS LOVE.DTA **CORR MATRIX** COR FACTORS" 

00251 /DA NI>20 NO«791 MA=KH 

00252 /MATRIX>IN(COR) 

00253 /MO NX=20 NK=2 PH=SY TD=SY 

00254 /LK 

00255 /MASCULUN FEMININE 



00256 /PA LX 

00257 /I 0 

00258 /I 0 

00259 /I 0 

00260 /I 0 

00261 /I 0 

00262 /I 0 

00263 /I 0 

00264 /I 0 

00265 /I 0 

00266 /I 0 

00267 /O 1 

00268 /O 1 

00269 /O 1 

00270 /O 1 

00271 /O 1 

00272 /O 1 

00273 /O 1 

00274 /O 1 

00275 /O 1 

00276 /O 1 

00277 /PA PH 

00278 /O 

00279 /I 0 



00280 /ST 0.6 LX(1,1) - LX(10,1) 

00281 /ST 0.6 LX(11,2) - LX(20,2) 

00282 /VA 1.0 PH(1,1) PH(2,2) 

00283 /OU SE MR MI FS SL>1 TM>1200 ND>3 

00284 
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